LOCAL AND GLOBAL STABILITY OF EQUILIBRIA FOR A CLASS 
OF CHEMICAL REACTION NETWORKS 



PETE DONNELL* AND MURAD BANAJpt 

Abstract. A class of chemical reaction networks is described with the property that each 
positive equilibrium is locally asymptotically stable relative to its stoichiometry class, an invariant 
subspace on which it lies. The reaction systems treated are characterised primarily by the existence 
of a certain factorisation of their stoichiometric matrix, and strong connectedness of an associated 
graph. Only very mild assumptions are made about the rates of reactions, and in particular, mass 
action kinetics are not assumed. In many cases, local asymptotic stability can be extended to global 
asymptotic stability of each positive equilibrium relative to its stoichiometry class. The results are 
proved via the construction of Liapunov functions whose existence follows from the fact that the 
reaction networks define monotone dynamical systems with increasing integrals. 
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1. Introduction. Systems of chemical reactions can give rise to dynamical sys- 
tems of various kinds (discrete or continuous time, discrete or continuous state, 
stochastic or deterministic, for example) and displaying a variety of behaviours [1]. 
Perhaps most widely studied are models whose evolution is naturally described by 
systems of ordinary differential equations, namely deterministic, continuous time, 
spatially homogeneous models where chemical concentrations take nonnegative real 
values. A broad question of interest is when models of some chemical reaction net- 
work (CRN) can be shown to allow, or forbid, certain behaviours for all reasonable 
choices of chemical reaction rates (kinetics). The attempt to make claims about the 
behaviour of CRNs which are to some degree independent of choices of kinetics is often 
termed chemical reaction network theory (CRNT). Since the pioneering work of 
Feinberg [5] and Horn and Jackson [3] , there has been considerable progress in various 
directions, including the discovery of structural features of networks associated with 
multistationarity, oscillation, and the persistence of solutions. 

Much, though not all, work in CRNT has focussed on reaction networks with 
mass action kinetics, but unknown rate constants, namely on particular polynomial 
differential equations with unknown parameters. Here we construct a class of CRNs 
which can be proved to have strong convergence properties with weaker assumptions 
on the kinetics. First note that the evolution of a CRN quite naturally takes place on 
certain invariant convex sets termed stoichiometry classes (to be defined below). The 
basic convergence properties of the networks we describe are: 

1. No more than one positive equilibrium on each stoichiometry class, and local 
asymptotic stability of each positive equilibrium on its stoichiometry class. 

2. Under additional assumptions, global asymptotic stability of a unique positive 
equilibrium on each nontrivial stoichiometry class. 

The precise meaning of these statements will be clarified below. The results 
will be proved using the theory of monotone dynamical systems [U [5]. There is a 
considerable intersection between this theory and the study of CRNs, reflecting the 
fact that CRNs fairly frequently give rise to order-preserving dynamical systems - see 
O [71 IHl [SI [ini [HI [m [I2 for example. The key geometrical insights for the results 
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presented here come from the results on monotone systems with increasing integrals 
in [13] , generalised in [Kl US] . 

2. Statement of the results. The local and global results summarised above 
will be stated precisely as Theorems 12.11 and 12.21 below after some terminology and 
notation are introduced. Define M"q to be the nonnegative orthant in K", i.e. 

= {x e M" : x^>0 for i = l,...,n}. 

Similarly 

M2;o = {x e M" : Xi<0 for i = l,...,n}. 

A vector in M"g will be referred to as nonnegative, while one in int(R"Q) will be 
termed positive. As chemical concentrations are necessarily nonnegative, R"q is the 
natural state-space for ODE models of systems of chemical reactions. 



We will be considering dynamical systems of the form: 

x = rv{x), (2.1) 

with the following assumptions: 

Al. X e r is an n X m matrix, u : [/ — 5- R™ is and is defined on some 

open neighbourhood U of IR>q. 

A2. The reaction rates v{x) satisfy conditions Kl, K2 and K3 listed in Ap- 
pendix |b] 

A3, r = AO, where 

i. A is an n X r matrix with each row containing exactly one nonzero entry, 
and no column of zeros. 

ii. is an r X TO matrix such that QijQkj £ for i ^ k and ker{Q"^) is one 
dimensional and includes a positive vector. 

A4. The DSR graph for the system at each x £ int(R"Q) is strongly connected. 

Remark on condition Al. x describes the concentrations of a set of n chemical 
species involved in to chemical reactions. The matrix F is termed the stoichiometric 
matrix of the system and Tij defines the net production/consumption of species i 
by reaction j. It is convenient, though not necessary, to assume that the vector 
field is defined on an open set containing the nonnegative orthant, in order to avoid 
technicalities when discussing its derivative. 

Remark on condition A2. This is a weak assumption on the kinetics. It 
can be crudely summarised via the statements: "reactions need all their reactants to 
proceed" and "provided all reactants are present, increasing a concentration speeds up 
a reaction" . The assumption has been discussed and illustrated previously in [TTl [13] 
for example. The assumption implies, amongst other things, that the nonnegative 
orthant is positively invariant (Lemma 10 in |13j). With condition Al, this guarantees 
that p.ip defines a local semiflow on R>g. 

Remark on condition A3. Condition A3 (i) implies that A has rank r, and 
Condition A3 (ii) implies that each column of Q contains exactly one negative entry 
and exactly one positive entry. The implications of condition A3 will be explored later 
during proof of the results. 

Remark on condition A4. The DSR graph associated with a CRN pTj is a 
signed, labelled, bipartite, multidigraph, with relationships to other well-known ob- 
jects such as Petri nets. It is strongly connected if there is a (directed) path from 
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each vertex to each other vertex. We need only the following reduced construction 
here. Given an n x m matrix A and an to x n matrix B, the (reduced) DSR graph 
Ga.b is defined as follows: it is a bipartite digraph on n + m vertices, Si, . . . , Sn and 
. . . , Rm with arc RjSi if, and only if, Aij ^ and arc SiRj if, and only if, Bji 7^ 0. 
If desired, the arcs may be given the signs of the associated entries in A and B. For 
System (|2.ip . at each x G R"q we can define G{x) = Gr^-Dv(x)- Under condition A2, 
G{x) is constant on int(M"g), and in fact on each elementary face of R>q (defined 
below) . 

Notation and terminology. Given any n x k matrix A, and x,y E M" wc will 
write X y for x — y E 1mA. Clearly is an equivalence relation on R". Given 
any such matrix A and any x G R>o, define 

Ca,. = {x + im{A)) n = {yeRlo -y 

In the study of chemical reactions, Cr,x is termed the stoichiometry class of x (also 
known as the "stoichiometric compatibility class" of x); for system (|2.ip satisfying 
assumptions A1-A4, Ca,x will be termed the A-class of x. Since A3 implies that 
Im(r) C Im(A), stoichiometry classes are subsets of A-classes, and clearly both are 
(forward) invariant under the local semifiow (p. Stoichiometry classes or A-classes 
intersecting int(R"Q) will be termed nontrivial. 

The first result of this paper is the following local claim: 

Theorem 2.1. Suppose that System satisfies assumptions Al, A2, A3 and 

A4-. Each equilibrium e e int(M"Q) is the unique equilibrium on its stoichiom,etry class 
Cr.e o,nd is locally asymptotically stable relative to Cr,e- 

This will be proved via construction of a Liapunov function on a neighbourhood 
in Cr,e of any positive equilibrium e. In order to extend this result to a global re- 
sult, we need conditions to ensure that each nontrivial stoichiometry class contains 
an equilibrium, the Liapunov function exists on the whole relative interior of each 
nontrivial stoichiometry class, and moreover that trajectories cannot approach the 
relative boundary of a nontrivial stoichiometry class. To make clear these notions we 
need to introduce some additional ideas. 

Given an n x r matrix A, define the closed, convex cone K{K) C R" associated 
with A as: 

K{K) ^{Ky.ye R^>o} . 
A local semifiow on IR"q is persistent if 

X e int(R|o) ^ ^{x) n 9R|o = 

where u){x) is the w-limit set of x. 

Let 5 C {1, . . . , n} be a subset of {!,..., n}, and 5^ = {1, . . . , n}\S. Define 

Fs = {a; G R" : x,>0,iClS and = 0, i ^ S}. 

Fs will be referred to as an elementary face of R"q. (Elementary faces are the rel- 
ative interiors of the closed faces of R>q.) An elementary face other than int(R"g) or 
{0} will be termed nontrivial. An elementary face Fs is repelling if, at each x G Fs 
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there exists « G S*"^ such that Xi > 0. Quite generally a repelling face of R"q can 
contain no w-limit points of a local semiflow on M"q (Lemma 11 in [13 ). 

To the list of conditions A1-A4 we add two further conditions: 

A5. K{A)nRlQ = {0}. 

A6. i. All reactions are reversible, or 

ii. Every nontrivial elementary face of R"q which intersects a nontrivial 
stoichiometry class of the system is repelling. 

Remark on condition A5. It will be shown later that condition A5 guarantees 
that the Liapunov function constructed via conditions Al~A4 extends to the entire 
relative interior of each nontrivial stoichiometry class. 

Remark on condition A6. It will be shown in Lemma 14.51 that conditions 
A1-A3 and A6(i) imply condition A6(ii) which, by the remarks above, implies the 
following: given any x £ int(R"Q) and y G Cr,x, then U!{y) fl 9R"g = 0, namely 
Cr,x n 9M"q contains no limit points of the local semiflow 0. Note that this is a 
stronger conclusion than persistence of 01^^ . 

We have the following global result: 

Theorem 2.2. Suppose that System \2. 1\) satisfies conditions Al, A 2, A3, A4, 
A5 and A6. Then each nontrivial stoichiometry class contains exactly one equilibrium, 
which lies in int(M"Q), and is globally asymptotically stable relative to its stoichiom- 
etry class. 

In physical terms chemical reaction networks satisfying conditions A1-A6 have 
very simple behaviour: different initial conditions on the same nontrivial stoichiometry 
class converge to the same positive equilibrium. 

3. Examples. The proofs of Theorems 12.11 and 12.21 are somewhat involved, so 
some examples of their application are provided first. 

3.1. Example 1. Consider the system of three reactions involving four chemicals 
A^B + C, B^D, C + D^A (3.1) 




and the associated differential equation (|2.ip with assumptions Al, A2. The stoichio- 
metric matrix, F admits a factorisation T = AQ as follows: 

f-1 
1 -1 
1 

V 1 

Note that A and Q fulfil conditions A3. The DSR graph at each interior point (drawn 
under assumption A2) is shown in Figure 13.11 and is clearly strongly connected. So 
condition A4 is satisfied. 

K{A) consists of all vectors of the form [a,b, —a, c]'^ where a,b,c > 0, and so 
K{A) n K"q = {0}, i.e., condition A5 holds. Finally, all reactions are assumed to be 
reversible, so condition A6 holds. By Theorem 12 . 2 1 each nontrivial stoichiometry class 
contains exactly one equilibrium which attracts the whole stoichiometry class. 

We remark that this example also fulfils the conditions in [IT] , from which global 
stability follows from a rather different argument. This is not the case for the next 
example. 
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Fig. 3.1. The DSR graph for reaction system \3. 1\) at any interior point in IR"g. Edge-labels 
are omitted. Negative edges are represented as dashed lines, while positive edges are shown with 
bold lines: only connectedness is important for the results here, and so edge signs are not needed; 
however including these makes the relationship between the DSR graph and the network of reactions 
fS. 1}} clearer to see. Pairs of antiparallel arcs of the same sign are represented as single undirected 
edges. 



3.2. Example 2. 

chemicals: 

A^B + C, 



Consider the following system of 4 chemical reactions on 5 



B ^ D, 



C + D^A, 



C + E^A. 



(3.2) 



As before, make assumptions Al, A2 about the kinetics. The stoichiometric matrix, 
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It is again easy to confirm that condition A3 is satisfied, and that the DSR graph at 
each interior point under assumption A2 (Figure 13. 2|) is strongly connected, so that 
condition A4 is satisfied. 




Fig. 3.2. The DSR graph for reaction svstem \3.2\ at each interior point o/R"q. Conventions 
are as in Fiaure VS. 1\ 

K{A) consists of all vectors of the form [a,b, — a, c, d]'^ where a,b,c,d > 0, and 
so K{A) n ]R"o = {0}, i.e., assumption A5 holds. Condition A6 again holds auto- 
matically as reactions have been assumed to be reversible. By Theorem 12.21 each 
nontrivial stoichiometry class contains exactly one equilibrium which attracts the 
whole stoichiometry class. 

To illustrate the case where reactions are not all reversible, we consider one further 
example. 

3.3. Example 3. We consider an example sometimes termed an "enzymatic 
futile cycle" 12. and whose biological importance is discussed in some detail in [T8] . 
Global stability in this example can be proven in many ways ( [12. for example), but 
is also an immediate consequence of the results in this paper. 
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The enzymatic futile cycle: 



Si+E ^ ESi S2 + E 

S2 + F ^ FS2 ^ Si+F 

has stoichiometric matrix which factorises as follows: 
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and DSR graph shown in Figure 13.31 



ESi- 



-> S2 



E 



Si ^ 



F • 



■FS2 



Fig. 3.3. The DSR graph for reaction svstem f3.3\ which is clearly strongly connected. Conven- 
tions are as in Fiaure V3.1\ Note that as a consequence of the irreversibility of some reactions, some 
edges in the DSR graph do not have an oppositely directed partner and so appear as directed. 



Conditions Al and A2 are assumed, while conditions A3-A5 are easy to confirm. 
On the other hand, Condition A6(ii) can be directly checked. The details involve the 
notions of siphons and associated faces of M>o, developed in Appendix [Cl and are 
presented in Appendix [P] 

4. Proofs of the results. 

4.1. Summary of the proofs. The proofs proceed roughly as follows. The key 
objects of interest are A-classes which are in general invariant r-dimensional objects in 
M". Conditions A1-A3 will imply that on each A-class is monotone with respect to 
the ordering defined by A'(A). Additionally, if condition A4 holds, then (j) is strongly 
monotone in the relative interior of each nontrivial A-class (in a sense to be made 
precise below). Condition A3 (ii) allows the construction of an increasing (linear) 
first integral on each A-class such that stoichiometry classes become the level sets 
of this function. One implication is that stoichiometry classes are unordered with 
respect to the ordering defined by AT (A). Via a construction closely related to those 
in [131 [T3], conditions A1-A4 allow the construction of a Liapunov function on each 
stoichiometry class in a neighbourhood of any positive equilibrium, strictly increasing 
along nontrivial trajectories. Adding condition A5 allows this Liapunov function to 
be extended to the entire relative interior of the stoichiometry class, while assumption 
A6 ensures that trajectories on nontrivial stoichiometry classes cannot have cj-limit 
sets intersecting 9E"g. While some of the arguments are minor modifications of 
arguments in jl5| . often they are simpler than the general arguments there, and so 
the presentation will be mostly self-contained. 
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4.2. Basics on cones and partial orders. A closed, convex and pointed cone 
K C R" (i.e., a closed, convex cone additionally satisfying K r\ {—K) = {0}) will be 
termed a CCP cone. If a CCP cone K has nonempty interior, then it will be termed 
proper [19| . Given an n x r matrix A with rank r, then K{h) is an r-dimensional 
CCP cone in generated by r extremal vectors (the columns of A). However, we 
can consider K{K) to be proper on any coset of Im A (the affine hull of K{K)) in the 
following sense: given z S R'' and c G R", the map z H> c + Az is a linear bijection 
from R*" to a coset of Im A, which maps the proper cone R>q to c + K{K). Via this 
identification, standard results on monotone dynamical systems can be lifted to cosets 
of Im A. 

The symbols <,>,<,>, 3> will refer to the standard partial ordering on R". 
For example, for a, 6 6 R", a <h means 6 — a £ R"q\{0}. When the ordering is that 
defined by some other cone if, the alternative symbols ^, ^, ^, > will be used. 
(Where these symbols are used, the cone in question will be clear from the context.) 
For example, given some CCP cone K C R" and a, & G R", a <b means h — a ^ K . 
Normally a < & means 6 — a G \\it{K)] here the meaning will be extended so that 
given any CCP cone K C R", a < & means 6 — a £ relint(i\r). 

Given a partial order < defined by a cone K, and some set X , if a;, y e X implies 
x 7^ y, we will say that "X is iiT-unordered" . 

4.3. A-classes are lattices. The lemmas to follow will show that the cones 
considered here define orderings which make each A-class a lattice. Given a, & G R*", 
define c = a A 6 by q = minjoi, hi}, and c = a V 6 by = maxja^, bi}. 

Lemma 4.1. Let K he an nx r matrix with rank r defining a cone K{A) C R". 
Let c G R" be an arbitrary vector. 

1. Consider x,y ^ W , z = x A y, x' = c + Ax, y' = c + Ay, and z' = c + Az. 
Then (i) z' < x' and z' < y' and (ii) Any 6' G c + Im(A) satisfying b' ^ x' 
and b' < y' , also satisfies b' ^ z' . In other words, z' is the infimum of x',y' 
in the order defined by K{A). 

2. Consider x,y G R*", z — x\/y,x' = c + Ax, y' = c + Ay, and z' — c + Az. 
Then (i) z' >: x' and z' y y' and (ii) Any b' E c + Im(A) satisfying b' >: x' 
and b' ^ y' , also satisfies b' ^ z' . In other words, z' is the supremum of x' ,y' 
in the order defined by K{A). 

Proof. Part 1 will proved. The proof of Part 2 is similar. Since z < x, x — z = 
p G R>o- So x' - z' = {c + Ax) - (c + Az) = Ap, i.e., z' < x' . Similarly z' <y' . 
Consider a vector 6' = c + A6 satisfying b' < x' and 6' ^ y'. i.e., 

x' - b' = c + Ax - {c + Ah) = Api and y' - b' = c + Ay - [c + Ab) = Ap2 , 

where pi,P2 G R>o- Multiplying each equation by any matrix A' such that A' A = / 
gives b < X and b < y, implying b < x A y = z, i.e. z — b = p^ G R>o. So 
^' - 6' = (c + Az) - (c + A6) = Ap3 G K{A), i.e. h' ^z'.U 

Notation. Given any A satisfying the assumptions of Lemma l4Tl and x,y € R" 
with X y, we write x a y for the infimum of x, y under the order defined by K{A). 
Similarly x v y refers to the supremum of x,y under the order defined by K{A). 
Clearly {x v y) {x a y) x y. 

The next lemma shows that additional assumptions on A ensure that the inter- 
section of a coset of Im(A) with any closed order interval under the standard ordering 
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(including, for example, M"q itself) is a lattice under the order defined by K{K). 

Lemma 4.2. Let A be annxr matrix of rank r, and such that each row contains 
no more than one positive entry and no more than one negative entry. Given c G M", 
x',y' G c + Im(A) and t G R"; 

1- If x\ y' > t, then x' a y' > t. 

2. Ifx',y' > t, then x' y y' > t. 

3. Ifx',y' < t, then x' a y' < t. 
4- Ifx',y' < t, then x' v y' < t. 

Proof. Choose x,y €W such that x' = c + Ax, y' = c + Ay. If row i of A contains 
a positive entry, then define p{i) by Aj p(j) > 0, and let — Aj p^j); otherwise let 
p(i) = 1 and ai — 0. Similarly, if row i of A contains a negative entry, then define 
q{i) by Aj < 0, and let f3i — — Aj otherwise let q{i) = 1 and /3i = 0. Note that 
for each i, ai,f3i > 0. We get 

{Ax)t = ajXp(i) - /3ja;g(i) and (A?/)^ = ai?/p(i) - Ay<j(i) ■ 

We prove Statements 1 and 2. Statements 3 and 4 follow analogously. 

1. Let z = X Ay so that, by Lemma [4.11 z' = c + Az = x' a y' . 

z[^ Ci + {Az)i = + min{xp(i) , ?/p(i) } - /^^ inm{xq(i),yq(i)} . 

So either 

z- = + aiXp(i) - PiXgi^i) = Ci + {Ax)i >ti, or 
z- = Q + a^yp^i) - Piyq(i) = q + {Ay)i > ti, or 

z[ = Q + a^Xp(i) - Piyq(^ > Q + {Ax)i > ti (because yq(^ < ^^(i)), or 
= Q + Q;jyp(i) - fiiXqi^if > Q + (Ay)i > (because < y,(j)). 

In every case, z- > ii, so z' > t. 

2. Let z = a; V y so that, by Lemma [4.11 z' ^ c + Az — x' v y' . 

z[ = Q + (Az)j = Q + max{a;p(i) , ?/p(i) } - max{a;q(i), ^^(i)} . 

So either 

z- = Q + aiXp(j) - = + (Ax)i > t^, or 

z- = Q + aj?;p(i) - = Q + (Ay)i > t;, or 

z- = Q + aia;p(i) - /3iy,(j) > c; + (Ay)^ > t., (because x^^^-^ > yp(^-,), or 
z- = Q + ai?/p(i) - ^iXq(i) >Ci + {Ax)i > ti (because yp(^o > Xp(i)). 



In every case, z^ > ii, so z' > t. □ 



Corollary 4.3. Let A be an n x r matrix with rank r, and such that each 
row contains no more than one positive entry and no more than one negative entry. 
Then, with the partial order defined by K{A), Ca,c is a lattice for each c G IR>o. 
Consequently, with assumption A3, each A-class of \2.1\) is a lattice. 

Proof. This follows from Statements 1 and 2 of Lemma 14.21 with t as the origin. □ 
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4.4. Each A-class has an infimum. We show that including condition A5 
ensures the existence of a unique minimal element on each A-class. 

Lemma 4.4. Consider any nxr matrix A of rank r, and such that K{K) HM^q = 
{0}. Then given any c G IR>q, Ca.c contains a greatest lower bound in the ordering 
determined by K{h). 

Proof. Consider any chain (i.e., totally ordered subset) C C C\^c- Take any 
decreasing sequence {xi) C C, i.e., xq xi >~ X2 >~ ■ ■ ■ ■ We want to show that 
(xi) is bounded (and hence bounded below). Since the sequence is arbitrary, this will 
imply that C is bounded below. Suppose, on the contrary, there exists a subsequence 
(xii^) with \xii^ I — )■ oo, implying that \xi^ — xo\ —5- oo. Consider the bounded sequence 
(xi^, — XQ)l\xi^\., and by passing to a subsequence if necessary, assume that it is 
convergent, i.e., (x^^ -a:o)/|xiJ y. But (xi^ -xo)/\xiJ = Xii^/\xiJ -a;o/|xiJ, and 
since Xi^/\xi^ \ is nonnegative and of magnitude 1, and xo/\xii^ \ — ?• 0, it follows that 
y is nonnegative and of magnitude 1. On the other hand, {xi^. — xo)/\xi^, \ £ —K{A), 
and since —K{A) is closed, y £ ~K{A). This implies that —K{A) includes a nonzero, 
nonnegative vector, contradicting the assumptions of the theorem. Thus every chain 
has a lower bound. By Zorn's Lemma, Ca,c contains a minimal element. Since Ca,c is 
a lattice, this minimal element is unique. □ 

4.5. With reversibility, oj-limit sets of nontrivial initial conditions lie in 

int(M"Q). The next claim is that if all reactions are reversible, no limit sets of System 
(j2.ip satisfying conditions A1-A3 on a nontrivial stoichiometry class can intersect 
9K"q. In other words, conditions A1-A3 and A6(i) imply condition A6(ii): 

Lemma 4.5. Consider System h2.1\) satisfying conditions A1-A3. Suppose ad- 
ditionally that all reactions are reversible, and let S be some nontrivial subset of 
{1, . . . If for some c ^ 0, Cr,c intersects Fs, then Fs is repelling. 

This claim involves a somewhat lengthy digression, and so is proved and illus- 
trated in Appendix [Cl Note that the assumption that chemical reactions are reversible 
means, mathematically, that the rate function v{x) fulfils conditions Kl and K3 in 
Appendix [B] 

4.6. An increasing first integral and unordered stoichiometry classes. 

In this subsection, make only the following assumptions on F, implied by assumption 
A3. Assume that F = AQ where 

Cl. A is an n X r matrix with rank r. 

C2. is an r X TO matrix, ker(0"^) is one dimensional, and there is a unit vector 
ye £ kcr(6-^) satisfying j/e ^ 0- 
Since A^ has rank r, it defines a surjective map from R" — s- W , and so we can choose 
and fix a vector pe such that A^pQ — ye- Note that pe G iiit(ii'(A)*), the dual cone 
to K{A), since for any z > 0, PqAz — y^z > 0. Define the linear scalar function 
H : M" ^ R by H{y) — p^y. The next lemma shows that H is increasing with 
respect to the order defined by K{A). 

Lemma 4.6. Consider € R" such that y > x. Then H{y) > H{x). 

Proof. Note that y = x -\- Az where z > 0. Then 

H{y) - H{x) ^ p%y - p%x = Pe{Az) = y^z > 0, 
where the last inequality follows because z > 0, and ye ^ 0. □ 
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Next, we prove that restricting attention to a A-class, stoichiometry classes are 
precisely the level sets of H . 

Lemma 4.7. For any x G Wi.^, Cr,x = {y^'^ x\ H{y) = H{x)}. 

Proof. From the definitions, Cr.x ^ Ca,x- Note that p^T = PqAQ — y^Q — 0. So 
given y x, and writing y = x + Ty' , 

H{y) = ply ^pI{x + Ty') = p^x - Hix) . 

SoCr,x C {y x\H{y) = Hix)}. 

On the other hand, consider any y ~^ x such that H{y) = H{x). Write y = 
X + Ay' and note that 

= H{y) ~ H{x) ^pl[y-x)^ plKy' - yly' . 

Since ye € ker(e^) and Im(e) = [ker(e^)]-L, y' e lm(e). So y' = Oy" for some y", 
i.e., y = x + AQy" = x + Ty" . Thus y x. So Cr.a: 2 {y x\ H{y) = H{x)}, 
proving the claim. □ 

Corollary 4.8. Each stoichiometry class of System V2. 1]) with assumption A3 
is K{ A) -unordered. 

Proof. Note that condition A3 implies conditions CI and C2. By Lemma 14.71 
X y implies H{x) = H{y). On the other hand, by Lemma 14.61 H{x) — H{y) 
implies x y, and it follows that Cr.y is A'(A)-unordered. □ 

Notation. Given the characterisation in Lemma 14.71 when we restrict attention 
to some A-class Ca_c, we can refer to the stoichiometry class in Ca.c on which H{-) 
takes the value h as Cjy ^, i.e. define Cjj ^ = {y G Ca,c | H{y) — h}. 

4.7. ii'-quasipositivity. Notation for matrices. Given any matrix M, we 
refer to the fcth column of M as Mk and the fcth row of M as . It is convenient to 
phrase results in terms of "qualitative classes" of matrices and related ideas. Given a 
matrix M, the matrix-sets Q{M), Qo{M) and Qi{M) are defined in Appendix Rl 

Terminology. Given a CCP cone K C R", and an n x n matrix J, we say that 
J is /iT-quasipositive if there exists a G M such that J + al : K K . If in fact 
there exists a e R such that J + al: K ^ K and J + al is i^- irreducible (namely, if 
F is a closed face of K and J + al: F ^ F, then either F — {0} or F = K), then we 
say that J is strictly iiT-quasipositive. 

The following lemma on ii'-quasipositivity of a special class of rank 1 matrices 
appears as Theorems 5.3 and 5.4 in [10 . 

Lemma 4.9. Consider a vector T, any vector V £ Qo(~r) and a CCP cone K 
with extremals {yi\. Define two conditions as follows: 

^. F = ryk for some k, and either (i) r > and yj G Qi(— L) for all j =/= k or 

(ii) r < and yj G Qi{T) for all j ^ k. 
B. There exist ri,r2 > and y, G Qi{T)\Qi{-T) , y^ G Qi(-F)\Qi(F) such 
that T = rij/i - r2yj. Moreover, yk G (Qi(F) n Qi(-F)) for k ^ {i,j}- 
If either F = 0, or A or B holds then TV'^ is K -quasipositive. 

Proof. The case F = is trivial. For the remaining cases, see [TU]. Although 
cones were assumed in that reference to be proper, the proofs are straightforward 
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calculations which apply for any CCP cone. □ 
An immediate corollary is: 

Corollary 4.10. Consider a matrix T, any matrix V G Qoi—^'^) and some 
CCP cone K . If each nonzero column of T satisfies either condition A or condition 
B of Lemma \4.9[ then TV is K -quasipositive. 

Proof. TV = X^i^iF*, and by Theorem 14.91 for each i there exists a.; such that 
T^V + a,I -.K^K. Clearly, defining a = Y,^a^,TV + al : K ^ K. U 

This leads to: 

Lemma 4.11. Let A he an n y. r matrix with rank r, and no more than one 
nonzero entry in each row. Let O he an r y. m matrix such that each column of O 
contains no more than one positive entry and no more than one negative entry. Let 
r = Ae,. V e Qo{-T^)- Then TV is K (A) -quasipositive. 

Proof. We need to show that for each i, Fj fulfils the conditions of Lemma [4.91 
The trivial case 8i = implies = 0. The reader can easily confirm that if 8i 
contains a single nonzero entry, then F^ satisfies condition A, and if Qi contains a 
positive entry and a negative entry, then F; satisfies condition B. □ 

Corollary 4.12. Consider System h2.1\) with assumptions A1-A3. At each 
X G R>o, the Jacohian matrix TDv{x) is K (A) -quasipositive. 

Proof. Assumption A2 implies that Dv{x) G Qo(~r^) at each x £ R>o- 

As- 
sumption A3 implies the assumptions of Lemma 14.111 which now gives the result. 
□ 

4.8. Strict iiT-quasipositivity. The aim in this section is to infer the following 
for System (|2.ip satisfying conditions A1-A3: at any point x where the DSR graph 
G{x) is strongly connected, the Jacobian J{x) is strictly i4r(A)-quasipositive (i.e., there 
exists a e M such that J{x) +a/: K{A) K{A) and J{x) -\-aI is if (A)-irreducible). 
The following is an immediate consequence of Theorem 1 in |20) . 

Lemma 4.13. Let K C M" he a CCP cone, A an nx m matrix, and B' an mx n 
matrix. Suppose that 1mA % spani^ for any nontrivial face F of K , and that AB is 
K -quasipositive for each B G Qq(B'). Then whenever the DSR graph Ga,b is strongly 
connected, AB is strictly K -quasipositive. 

We wish to apply this resuh with A = T, B' = -F^, and K = K{A). We have 
already seen that AB is i^T-quasipositive for all B G Qo(— T^), and so we need: 

Lemma 4.14. Consider an n x r matrix A with rank r and an r x m matrix 
Q with no row of zeros. Let T = AQ. Then ImF ^ span.F for any nontrivial face 
F C K{A). 

Proof. Since rank A = r, a vector z G ImA has a unique representation z = 
^a^Ai. Consider some nontrivial face F of K{A), and choose some k such that 
Afc ^ F. Since Q has no row of zeros, choose some i{k) such that Ok.iik) — a 0. 
Define y = Tei(^k) = ^^(^i{k) — crA^ + • • • . Since this representation of y is unique, 
clearly y ^ spani^ and so ImF ^ spani^. □ 

We can deduce that: 

Corollary 4.15. Consider System 12. 1]) with assumptions Al-AS. Assume that 
at some x G ]R"q, the DSR graph Gr.Du(x) strongly connected. Then the Jacohian 
matrix TDv{x) is strictly K (A) -quasipositive. 
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Proof. Condition A3(ii) implies that Q has no row of zeros, for if Qij — for some 
i and all j, and z is a positive vector in ker which exists by assumption, then {z, e^} 
are two linearly independent vectors in herQ^. So by Lemma 14.141 ImF ^ span_F 
for any nontrivial face F C K{A). Since, by Lemma |4TTT1 rDv{x) is ii'-quasipositive 
for all Dv{x) £ Qo(^r-^), Lemma [4.131 now implies that TDv{x) is strictly -f^(A)- 
quasipositive whenever G'r,_Dt,(a;) is strongly connected. □ 

It follows immediately that: 

Corollary 4.16. Consider System \2.1\) with assumptions A1-A4-. The Jaco- 
bian matrix TDv{x) is strictly K (A) -quasipositive at each x G int(R"Q). 

4.9. Monotonicity and strong monotonicity. We have shown in Corollar- 
ies 14.121 and 14.161 that conditions A1-A4 imply that TDv is i^T (A)-quasipositive on 
all of R"q and strictly i4r(A)-quasipositive on int(R"Q). The implications in terms of 
monotonicity and strong monotonicity of the local semiflow restricted to A-classes are 
discussed briefly. 

Notation. Given a, 6 £ M", define the closed segment [a,b] ~ {\a + (1 — X)b : 
A G [0, 1]}. Open and semi-open segments (a, &), (a, b] and [a, b) are similarly defined. 

The following is a consequence of results in [4] (see also |13j): 
Lemma 4.17. Consider a proper cone K C W , some open set U , and a 

vector field f : U ^ W with Jacohian matrix Df. Let X C U be some convex 
domain, positively invariant under the local flow (pjj defined by f , and let cj) be the 
induced local semiflow on X. Assume that Df is K -quasipositive in X. Consider 
some x,y d X with x <y. Then 4>t{x) ^ 0t(y) for each t > such that (j)t{x) , (j)t{y) 
exist. If X ^ y and there exists z G such that Df{z) is strictly K -quasipositive, 

then 4>t{x) -< 0t(y) for each t>0 such that (j)t{x) ^ <j)t{y) exist. 

Our interest is in A-classes. Recalling that A can be regarded as a bijection from 
W to Im(A), mapping R>g to Ki^K), Lemma 14.171 implies : 

Corollary 4.18. Consider System h2.1\) with assumptions A1-A4 and x,y G 
R"q with X <y. Then (j^tix) ^ 'f'tiy) for each t > Q such that <pt{x) , iptiy) exist. If 
X < y and at least one of x,y is in int(R"Q), then 4it{x) 0t(2/) (namely, (j)t{y) — 
4>t{x) G relint _fi'(A) j for each t> Q such that (j)t{x) , <j)t{y) exist. 

4.10. Structure of the equilibrium set. Define E C R"(, to be the equilib- 
rium set of (|2.ip . i.e. 

E^{x^ R'^o : Tv{x) = 0}. 

Lemma 4.19. Consider System i2. 1\) with assumptions A1-A4 and two distinct 
equilibria x,y with x y and at least one of x,y G int(R"g). Then either x y or 
y ^ X. Consequently no stoichiometry class with an equilibrium in int(R"Q) contains 
more than one equilibrium. 

Proof. Assume, by relabelling x and y if necessary, that x )/- y. 

(i) Suppose X < y, but x y. Since at least one of x or y lies in int(R"g), the 
line segment [x,y] certainly intersects int(R"Q). CoroUarv 14.181 then implies 
that for t > 0, a; = 4>tix) < <ptiy) = y, 3. contradiction. 

(ii) Now suppose that x and y are unordered. Then z = x Ay is different from x 
and y. By monotonicity (ptiz) ^ (t^tix) = x and ^((z) ^ 0t(y) — y- Since z = 
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X Ay, 4>t{z) :< z. But (f)t{z) e Cr.z by invariance ofCr.z, and Cr,z is unordered 
(Corollary S^, so z e E. If y e int(IRt^o) (resp. x G int(R?io)), applying the 
argument in part (i) to z and y (resp. z and x) gives a contradiction. 
Since stoichiometry classes are subsets of A-classes and are unordered (Corol- 
lary |4]8]), it follows immediately that no stoichiometry class with an equilibrium in 
int(K"Q) contains more than one equilibrium. □ 

Remark. Strict ordering of equilibria followed from a considerably more involved 
argument in [15j . Here the fact that A-classes are lattices makes the conclusion simple. 

Given x y, define 

P(x, y) = {{x + if (A)) U{x- K{A))) n Cr,y • 
Using the characterisation in Lemma 14.71 we can alternatively write 

P{x, y) = {{x + K{K)) U (x - K{K))) n {z a; : H{z) = H{y)}. 

Lemma 4.20. Assume that matrices A, and F satisfy Condition A3. Then 
P{x, y) is nonempty for any y x. 

Proof. If H{y) = H{x), then by Lemma |4.7[ x e Cr.j, and we are done. Assume 
that H{y) > H{x) (resp. H{y) < H{x)). By the lattice property of Ca,i^ (Corol- 
lary US]), z — X Y y C\^x n (a; + if (A)) (resp. z — x a y G Ca,k fl (x — K[A))), 
and by Lemma H{z) > H{y) > H{x) (resp. H{z) < H{y) < H{x)). Since 
{x + if (A)) n Ca,x (resp. {x ~ if (A)) n Ca,x) is convex, it includes [x,z]. By the 
intermediate value theorem there exists w £ [x,z] such that H{w) ~ H{y). By con- 
struction w G P{x,y). □ 

Lemma 4.21. Consider any CCP cone if C ]R"g and some vector p € int(if*). 
Civen any point x £ R" and any constant t > 0, the set {x + K)r\{y : p^y — p^ x-\-t^ 
is hounded. Similarly, [x — if) H {y : p^y = p^ x — t} is bounded. 

Proof. The first statement will be proved; the proof of the second is similar. With 
fixed X and p, define 

R — inf p^y. 

yGK,\y\ = l 

Since p^ {y — x) > (as p G if* and y — xE if \{0}), i? > as the infimum of a positive 
function on a compact set. Consider any sequence y„ in {x + K)r\{y : p^y = p^x+t}. 
We then have i = p-^(y„ —x) > R\yn — x\, i.e., \yn — x\ <t/R, and so (y„) is bounded. □ 

Lemma 4.22. Assume that matrices A, O and T satisfy Condition A3. P{x,y) 
is a nonempty compact, convex set for any y x. 

Proof. It has been shown in Lemma 14.201 that P{x, y) is nonempty. Either x = y, 
in which case P{x, y) = {x} which is trivially compact and convex; or exactly one of 
{x + if (A)) n Cr,y or (x — if (A)) n Cr,y is nonempty (since stoichiometry classes are 
unordered by Corollarv l4.8p . For definiteness assume that P{x,y) = {x + K{A))nCr,y 
is nonempty (the other case is similar). P{x,y) is thus closed and convex as the in- 
tersection of closed, convex sets. Applying Lemma [4.211 with if — if (A) and p — pq, 
P{x, y) is bounded. □ 

Lemma 4.23. Consider System \2.1]) satisfying conditions A1-A3. Choose some 
e £ E. Then each stoichiometry class in C\_e contains an equilibrium. 
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Proof. By Lemma 14.221 for arbitrary x e, the intersections P(e, x) are 
nonempty, compact, convex sets. Moreover these sets are forward invariant under 
((e + K{A)) U (e — K{A))) is forward invariant by monotonicity of cj) (Corollary 14. ISp 
and since the stoichiometry class Cr,x is forward invariant, P{e, x) is forward invariant 
as the intersection of forward invariant sets. By the Brouwer fixed point theorem f |21) 
for example) P{e,x) contains an equilibrium. Since x was arbitrary, each stoichiom- 
etry class in Ca.g contains an equilibrium. □ 

Lemma 4.24. Consider System \2.1\) satisfying conditions Al-A^ and some e G 
E n int(M"Q). There exists an equilibrium eo S (e — relint(iir (A))) fl int(M"g). Given 
any such equilibrium cq there exists a homeomorphism 

^ : [H{eo),H[e)]^^{[H{eo),H{e)]) C E 

such that (i) H{tp{h)) = h, (ii) V-Ci-ff (eo), i?(e)]) C int(R'^o) and (Hi) hi < h2 ^ 
ip{hi) <il:{h2). 

Proof. Certainly, e is a not a minimal element of CA,e since it lies in int(IR"Q), 
and moreover e lies in relint(Cr,e)- Define R as in the proof of Lemma I4.21[ fixing 
K ~ K{A) and p — pQ, namely: 

R = inf PpiV- 

Let < (5 be the minimum distance from e to (9R"q and choose some positive e < 5R. 
Then for any x < e such that H{e) — H{x) < e, P(e, x) C int(R"Q): this follows since 
\y - e\ < (i?(e) - H{x))/R < S ior y & P{e,x) (see the proof ~of Lemma li?^ . By 
Lemma [4.221 P(e, x) is a nonempty, compact convex set, and consequently contains an 
equilibrium eH(x)- Since eH(x) S (e — relint(iir(A))) n int(R"Q), it is, by Lemma [4. 191 
the unique equilibrium in Cr,x- Defining eo = eH{e)^e, for h G [H{eo), H{e)], the map 
ip : h Ch is thus well defined, has image in int(R"Q) and is clearly a bijection onto 
its image. is continuous as it is simply the restriction of H to the image of ip. 

We next show that ip is continuous (see also the proof of Lemma 5.12 in [TS]) 
Consider any h G [H{eo), H{e)], a sequence of values hi C [H{eo),H{e)] with hi — >■ h 
and the corresponding equilibria = ip{hi). Since all Ci lie in the order interval 
[[eo,e]] = {?/ G K" : cq ^ y ^ e} C int(]R"o) which is easily seen to be bounded, e 
contains no unbounded subsequences. Consider any convergent subsequence of (e^) 
say ei^ e. By closure of E and [[eo,e]], e e E Ci [[eo,e]], and by continuity of 
H, H{e) = h. Since V'(^) is the unique equilibrium satisfying these requirements 
e = ip{h). Thus ip is continuous. 

Finally, that H{tp{h)) = ft, is immediate from the definition, and that /ii < /i2 
tp{hi) -€ ip{h2) is immediate from Lemma [4.191 since lm{tp) C int(M"Q). □ 

Remark. We could equally prove the existence of eg G (e + relint(i4r(A))) fl 
int(E"Q) and a homeomorphism -0 : [i?(e), i?(eo)] — J> ip{[H{e), H{eQ)]) C E. 

4.11. Local asymptotic stability of equilibria. We are now in a position to 
prove the local stability of all positive equilibria for System (|2.ip satisfying conditions 
A1-A4. 

Proof of Theorem \2.1\ Consider any positive equilibrium e. That e is the unique 
equilibrium on Cr.e follows from Lemma 14.191 Choose and fix some equilibrium eg G 
(e-relint(i^(A)))nint(M"p) as in the proof of Lemma 14. 241 Via Lemma define a 
strictly increasing homeomorphism t/j : [iJ(eo), -ff (e)] — > ip{[H [ei^) , H [e)]) = -Eeo,e ^ E 
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such that ip{H{ea)) — eg and tp(H{e)) — e. Note that £'eo,e Q int(K>Q) and is strictly 
ordered. Define U ^ (cq + K{A)) n Cr.e- For any point x e U\{e}, cq e x - K{A). 
On the other hand e ^ x ~ K{A). Since -Eeo,e is homeomorphic to a line segment (and 
hence connected), Eeg.e must intersect x — relbdiir(A). Moreover this intersection 
is unique, otherwise some pair of distinct equilibria in i?eo,e must fail to be strictly 
ordered. It is also clear that (e — relbd-fir(A))ni?e(,^e = {e}. Thus for all x ^ U, define 
Q{x) = {x — relbdi^(A)) n i?eo,e, and L{x) = H{Q{x)). We make three claims (the 
reader may wish to compare Lemmas 5.16, 5.17 and 5.18 in |15)): 

1. L(e) > L{x) for x G /7\{e}; 

2. Q, and hence L, are continuous on U ; 

3. L increases strictly along nontrivial orbits. 

The first statement is immediate since H{e) > H{z) for z G i?e(,^e\{e}. Since both 
K(A) and £'eo,e are closed sets, with intersection at a unique point, it is not hard to see 
that for (xi) C U with Xi x, Q{xi) — > Q{x) and the second claim follows. Finally, 
given X G U\{e}, by strong monotonicity (Corollarv l4.18p . (f>t{x) > <j)t{Q{x)) = Q{x) 
for t > 0; if z ^ Q{x) then z -C 0t(a;), i.e., z ^ Q{(l)t{x)). So Q{(l)t{x)) >- Q{x), and 
thus L{(j)t{x)) > L{x), proving the third claim. 

Thus L serves as a Liapunov function for and, by standard arguments, e 

is locally asymptotically stable relative to Cr,e- In particular, each x G relintC/ is 
attracted to e. □ 

4.12. Global asymptotic stability of equilibria. For a global result we need 
to show that the previous local constructions can be extended. 

Lemma 4.25. Consider System V2. 1]) with conditions A1-A4 defining a local 
semiflow <j). Consider some c G int(R"Q). Suppose Cf^,c contains an infimum z ^ C\^c, 
and that given any y G Ca,c H int(IR"g), has no Lu-limit points on Cr^y H i9]R"q. 
Then, for each h G [H{z),H{c)], the stoichiometry class contains exactly one 
equilibrium. For h G {H{z),H{c)] this equilibrium is in int(R"Q). 

Proof. Since stoichiometry classes are unordered (Corollarv 14. 8p . and z ^ Ca,c, 
{z} must be an entire stoichiometry class, and consequently (since stoichiometry 
classes are invariant) z is an equilibrium. Since Ca,c contains an equilibrium, by 
Lemma 14.231 each stoichiometry class in C\_c contains an equilibrium. Since Cr.c is a 
nontrivial stoichiometry class, by the assumption on cj-limit sets of </> all equilibria in 
Cr,c lie in int(R"Q). By Lemma |4.19[ there is in fact a unique equilibrium in Cr.c- In 
summary, Cr,c contains a unique equilibrium and this equilibrium lies in int(R"Q). 

Choose any h G {H{z), H{c)]. By continuity of H, there exists x G {z,c] such 
that H{x) — h, i.e. x € C'^ ^. By basic properties of convex sets, x G int(R"Q), and 
so Cr,x is nontrivial. Applying to Cr.x the argument applied to Cr.c, Cr,x contains 
exactly one equilibrium, and this equilibrium is in int(R"Q). □ 

Lemma 4.26. Consider any c G int(R"Q) and the nontrivial A-class Ca,c with 
the assumptions of Lemma \4-25\ There exists a strictly increasing homeomorphism 
ip : [H{z),H{c)] iJj{[H{z),H{c)]) C E and such that (i) H{i;{h)) = h, (ii) hi < 
h2 =^ ip{hi) < i){h2). 

Proof. By Lemma [4.251 for each h G [H{z)^ H{c)\, contains a unique equilib- 
rium e/i and for h G {H{z), H{c)], Ch G int(R"Q). Defining ij; by V'(^) — it is clear 
that V' is a bijection. Continuity of ip and -0"^ now follow as in Lemma 14.241 That 
hi < h2 ^ i^{hi) < '^'{^2) follows from Lemma 14.191 □ 



We are now finally ready to prove the main global convergence result. 
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Proof of Theorem \2.2\ Consider any c G int(M"Q), the nontrivial stoichiometry 
class Cr.c and the nontrivial A-class Ca^c- By Lemma [4.4[ assumptions A3 and A5 
imply that there exists z G Ca,c with z ^ C\^c- By condition A6, for y e int(M"Q), 
(j> has no w-limit points on Cr,y n 9]R"q. Thus Lemma [4.251 applies: for each h G 
[H{z), H{c)], contains exactly one equilibrium eu and for h G (i? (z), iJ(c)], Ch G 
int(]R"Q). Let e = eH{c) be the equilibrium in Cr,c- 

By Lemma[121 the map V : [H{z),H{c)] tp{[H{z), H{c)]) = E^^^ Q E defined 
by ipih) — Ch is a, strictly increasing homeomorphism. We now follow the proof of 
Theorem [O with eo = z. Since z ^ Ca,c, P{z,c) = {z + K{A)) n Cr.c = Cr,c- 
One immediate consequence is that Cr,c is bounded by Lemma 14.221 Thus given any 
y G Cr,c, "^(y) (the w-limit set of y) exists. 

For y e Cr,c we define Q{y) — [y ^ relbdX) fl E^^c and construct the Liapunov 
function L(-) = H{Q{-)) defined on all of Cr,c- That Q is well defined, L{e) > L{y) 
for y G Cr.c\{e}, Q, and hence L, are continuous on Cr,c follow as in the proof of 
Theorem 12.11 Note that L is strictly increasing, namely L{(f)t{y)) > L{y) for t > 0, 
if y G (Cr,c\{e}) n int(R5o), but if y G Cr,c n dW^^ we can only claim that L is 
nondecreasing, namely L{4>t{y)) > L{y) for i > 0: it may occur that Q{y) = z, in 
which case the entire segment [z, y\ may lie in 9M" and wc only have 4>t{y) ^ 4't{z) = z. 
Thus u}{y) lies in {e} U (Cr.c H 9M>q). But by condition A6, (j) has no w-limit points 
on 9K>0' ^^'^ ^iv) — {^}- Thus all initial conditions in Cr,c are attracted to e. □ 

5. Conclusions. A class of CRNs has been described with strong convergence 
properties under only weak kinetic assumptions. The networks in this class are defined 
primarily by the existence of a certain factorisation of their Jacobian matrices, and 
strong connectedness of their DSR graphs. Roughly speaking, the convergence prop- 
erties of these CRNs spring from the combination of monotonicity of the associated 
dynamical systems and the existence of integrals of motion. 

For the class discussed here, monotonicity is with respect to an order defined by a 
cone with linearly independent extremal vectors. An interesting question is whether 
it is possible to extend the theory to more general cones thus progressing with the 
program of identifying CRNs with simple behaviour: while several of the proofs here 
were simplified by the fact that A-classes were lattices, the results in [15l US] suggest 
that this may not be crucial to the geometric argument. 

Finally, it was observed that Example 1 fell into a class which can also be proved 
to be monotone in "reaction coordinates" An interesting theme for future work 
is to work towards a synthesis of the approaches to monotonicity of CRNs in normal 
("species") coordinates and reaction coordinates. 
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Appendix A. Qualitative classes of matrices. A real matrix M determines 
the qualitative class Q{M) consisting of all matrices with the same sign pattern 
as M . Explicitly, X G Q{M) if and only if i) X has the same dimensions as A/, ii) 
[Mij > 0) {Xij > 0), in) {Mij < 0) (Xy < 0) and iv) (Afy = 0) ^ (X^j = 0). 

The closure of Q(M) is denoted by QoiM), namely X G QaiM) if and only if i) X 
has the same dimensions as M, ii) (Mij > 0) {Xij > 0), iii) {Mij < 0) => (Xij < 0) 
and iv) (Afy = 0) (Xy ^ 0). 
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Qi{M) is defined by deleting iv) from the defining properties of Qq{M). 



Appendix B. Kinetic assumptions. Closely following [13] the following natu- 
ral assumptions are made about the kinetics, namely about the function v{x) in (j2.ip . 
Ij^i is the set of indices of chemicals occurring on the left of reaction j and Ij^r is the 
set of indices of the chemicals occurring on the right of reaction j. 

(Kl) Tij{dvj / dxi) < 0, and if Tij = then dvj/dxi — 0. More briefly, in the 

e Qoi-T^) at each x. This 
for example, and is satisfied by 



dvi (x) 



notation of the previous appendix, 

condition has been discussed before in [22l |£ 
all reasonable kinetics provided no chemical occurs on both sides of a reaction. 
(K2) If reaction j is irreversible then 

(i) Vj > with Vj = if and only if = for some i S Tjj . 

(ii) If Xi > for all i £ Ijj, then dvj/dxi > for each i e Xjj. 
(K3) If reaction j is reversible then 

(i) If = for some i e X,_; (resp. for some i e 2j.r) then Vj < (resp. 
v,>0). 

(ii) If Xi = for some i G Ijj (resp. for some i G 2j.r), then vj < (resp. 
Vj > 0) if and only if Xi' > for each i' G Ij^r (resp. for each i' G 

(iii) If Xi > for all i G Ij^i (resp. for all i G Ij,r) then, for each i G X,,; 



(resp. i G Ij,r), dvj{x)/dxi > (resp. dvj{x)/dxi < 0). 
One consequence of K3(iii) is that for a reversible system of reactions 
Q(-rT) at each x G int(M?!.o). 



dxj 



Appendix C. Persistence in the reversible case. Given nonempty S C 
{1, . . . ,n}, define as the n x n projection matrix which maps the nonnegative 
orthant to Fs, the closure of the elementary face Fs Q M>q (i.e., Pfj = 1 if i = j G 5 
and Pfj — otherwise). A mixed-column matrix is a matrix such that each nonzero 
column contains both a positive and a negative entry. We note that as shown in [13j 
condition A2 implies that: 

1. An elementary face Fs of M"q is cither repelling or the vector field is every- 
where tangent to Fs- In the latter case the set 5'^ is termed a siphon [231124] . 
The face Fgc associated with a siphon S will be termed a siphon face. 

2. For a system of reversible reactions with no chemical appearing on both sides 
of any reaction, Fs fails to be repelling if and only if (/ — P^)T is a mixed- 
column matrix. In intuitive terms, if the concentration of chemical i from 
one side of reversible reaction j is zero on Fs, while concentrations of all 
chemicals on the other side of reaction j are nonzero on Fs, then Fs must be 
repelling since chemical i is being produced by reaction j. This occurrence 
manifests as (/ — P^)r failing to be a mixed-column matrix. 

Lemma C.l. Let S be a nonempty proper subset of {1, ... ,n}. Given an n x m 
matrix T, suppose there exists w satisfying V"'" [I — P^)w = and (/ — P^)w > 0. 
Then for each c G int(K"o)' ^r.c does not intersect Fs. 

Proof. Cr,c intersects Fs iff there exist y G K™, z G M>q which solve the equation 

c + Ty = P^z. 

Left- multiplying both sides of the above equation by uF {I — P^) gives 

w^(/ - P'^)c + uF{I - P^)Ty ^ uF{I - P^)c = vF{I - P^)P^z = 
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Since c » and (/ — P^)w > this is a contradiction, and so Cr,c H Fs = 0. □ 



The above lemma tells us that if there exists w such that (/ — P^)w > and 
r^(/ - P^)w = the n no nontrivial stoichiometry class intersects Fs- 



Lemma C.2. Suppose that A and are matrices satisfying the conditions in 
A3, and Q is a projection which takes ]R"q to some nontrivial (closed) face o/]R"g. 
If J = QAQ is a mixed-column matrix, then there exists w satisfying Qw > and 
e^A^Qw = 0. 

Proof. Suppose some column of A^Q is a nonzero multiple of a unit vector Ck and 
another is an oppositely signed multiple of Cfc. In particular, suppose {A'^Q)i = tiCk, 
(A'^Q)j ~ —t2ek, where ti,t2 > 0. Clearly Qii = ii and Qij = ej (otherwise 
(A^(5)i = 0, resp. (A^Q)^ = 0), and the vector w = t2ei + ticj satisfies Qw = w > 
and A^Qw = t2tiek — ^1^26^ = 0. 

Suppose K{A^Q) D Rt;^ or K{A'^Q) D W^„. If K{A'^Q) D MC^g, then for some 
w; > 0, A^Qw = ye, where <§; j/e G ker0^ and so Q^A^Qw — 0. Similarly if 
K{A^Q) 3 M<o, then for some w > 0, A^Qw = -ye, and so Q'^A^Qw = 0. 

Thus if there exists no w satisfying Qw > 0, Q^A'^Qw = 0, then each of the 
following must be true: 

1. Given a column of A^Q of the form tiSj where ti ^ 0, there is no column of 
A^Q of the form t2ej with tit2 < 0. 

2. K{A^Q) 2 KS:o and K{A^Q) ^ K!co- 

These together imply that each column of QA has fixed sign (i.e. it is either zero, or 
has some nonzero entries all of which are positive, or has some nonzero entries all of 
which are negative), and there exist two columns with different sign. More precisely, 
defining ^0, S+,S- C {1, . . . , r} by z G 5o O (QA), = 0, i G ^+ ^ {QA), G M|o\{0}. 
and i e S- ^ {QA)i G K<o\{0}, then {S'o, 5+, S*-} is a partition of {l,...,r} and 
at least two of So,S+,S- are nonempty. By the definition of Q there exists some 
column of Q, say column fc, containing nonzero entries Qik and Qjk such that i,j 
are not both in the same member of the partition {Sq, S+,S-} (i.e., such that {QA)i 
and {QA)j are of different sign in the sense defined above). If this were not the case, 
then Q could be written in block-diagonal form, in which case it is easy to show 
that ker O-^ contains a nonnegative vector which is not a multiple of ye . But then 
((5A6)fc = Qik{QA)i + Qjk{QA)j is a nonzero vector all of whose nonzero entries are 
of the same sign, and hence QAQ fails to be mixed-column matrix. □ 

Proof of Lemma \4^ If Fs fails to be repelling and the reactions are reversible, 
then (/ — P^)T is a mixed-column matrix. Applying Lemma [0.21 with Q = I — P^ , 
there exists w satisfying (/ — P^)w > and r-^(/ — P^)w = 0. But then, by 
Lemma IC. 11 Cr,c does not intersect Fs for any c 3> 0. In other words no nontrivial 
stoichiometry class intersects Fs- □ 



We illustrate the preceding lemmas with an example. Consider the matrices 



A: 



/ 1 \ 

1 

1 

-10 

V 1 y 



e 



1 





1 


1 


1 








-2 


-2 
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F = Ae = 



/ 


1 





1 \ 




-1 


1 







1 





1 




-1 





-1 


V 





-2 


-2 ; 



Consider the pairs 



Si 


= {2,3,5} 


Wl = 


(1,0,0,1,0)^ 


S2 


= {1,2,5} 


W2 = 


(0,0,1,1,0)^ 


S3 


= {2,5} 


W3 = 


(1,0,1,2,0)^ 


Si 


= {3,4} 


W4 — 


(2,2,0,0,1)^ 


S5 


= {1,4} 


W5 = 


(0,2,2,0,1)^ 


Se 


= {4} 


Wq = 


(1,2,1,0,1^ 


S7 


= {3} 


W7 = 


(3,2,0,1,1)^ 


Ss 


= {1} 


ws = 


(0,2,3,1,1)^ 



For each Si C {1,2,3,4,5} listed above, (/ — P^^ jT is a mixed-column matrix, and 
these are the only proper nonempty subsets of {1, 2, 3, 4, 5} with this property. Con- 
sider for example the set Si = {2, 3, 5}, giving 



/ 1 











\ 


/ 


1 





1 \ 


( 


1 





1 \ 



















-1 


1 

































1 





1 






















1 







-1 





-1 




-1 





-1 


V 











/ 


V 





-2 


-2 J 


\ 








/ 



which is clearly a mixed-column matrix. Further, in each case Wi satisfies (/ — 
P^')wi > and r^(/ - P^^)wi = 0. Thus, as predicted by Lemma [Q!2l when- 
ever (/ — P^)r is a mixed-column matrix, there exists a corresponding vector r, 
satisfying (/ — P^)r > and r-^(/ — P^)r = 0. Note that in the first three cases 
(/ - P'^'jw, G kcr A^, while in the rest A^(/ - P''^')w, = (2, 2, 1)^ e kerO^. 

Appendix D. Checking condition A6(ii) for Example 3. 

For the class of kinetics treated here, as described in [23], siphons can be identified 
with sets of chemicals for which each "input reaction" is also an "output reaction" ; in 
other words 5 is a siphon if and only if for each i G S every reaction able to produce 
species i requires some chemical j € S. Similar to the process in |23j using Petri 
net graphs, siphons can readily be identified from examination of the DSR graph, 
provided that arcs are signed. 

A siphon is minimal if it does not contain (strictly) any other siphon. Clearly 
every siphon face lies in the closure of at least one siphon face associated with a 
minimal siphon: so if we can identify each minimal siphon S and prove that nontrivial 
stoichionietry classes do not intersect Fs" , this ensures that nontrivial stoichiometry 
classes intersect no siphon faces at all. 

Recall that for Example 3 the stoichiometric matrix takes the form: 

1 \ 

1 
1 
1 -1 
-1 1 
1 -1 / 

corresponding to the substrate ordering [Si, E, ESi, S2, F, FS2]- In this case there 
are exactly three minimal siphons: Ei — {2, 3} corresponding to substrates {E, ESi}; 
S2 = {5, 6} corresponding to substrates {F, FS2}; and S3 — {1, 3, 4, 6} corresponding 
to substrates {5*1, ESi, S2, FS2}- To proceed, note that some nontrivial stoichiometry 



/ -1 

-1 

1 
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class intersects the closed face F^c if and only if there exists c G int(M>g) and z G M>g 
such that 

c + Ty = P^>z. 

However if such c G int(]R>Q) and z G ]R>q exist, then defining v\ — (0, 1, 1, 0, 0, 0)"^, 
V2 — (0, 0, 0, 0, 1, 1)"^ and i>3 = (1, 0, 1, 1, 0, 1)"^, we obtain the following contradictions: 

< vjc = vjc + ry = P^' z = (i = 1, 2, 3). 

These contradictions confirm that the equations c + Fy = P^* z cannot be satisfied 
for any c ^ and z > 0, and thus no nontrivial stoichiometry class intersects any of 
the three closed faces Pec. 
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